Se NMR measurements of the n — d exchange field in the organic conductor 

A-(BETS)2FeCl4 



Guoqing Wu,^ W. G. Clark/ S. E. Brown/ J. S. Brooks, ^ A. Kobayashi,^ and H. Kobayashi * 
^Department of Physics and Astronomy, UCLA, Los Angeles, California 90095, USA 
^National High Field Laboratory and Florida State University, Tallahassee, Florida 32306, USA 
"^Research Center of Spectrochememistry, University of Tokyo, Japan and 
Institute of Molecular Science, Okazaki, Japan 
(Dated: February 1, 2008) 

^'^Se-NMR spectrum and frequency shift measurements in the paramagnetic metal (PM) and 
antiferromagnetic insulating (AFI) phases are reported for a small single crystal of the organic 
conductor A— (BETS)2FeCl4 as a function of temperature (T) and field alignment for an applied 
magnetic field -Bo = 9 T. The results show that in the low T limit, where the localized Fe'^+ spins 
{Sd = 5/2) are almost fully polarized, the conduction electrons (Se yr-electrons, spin = 1/2) in 
the BETS molecules experience an exchange field (B.^£j) from the Fe'^^ spins with a value of — 32.7 
± 1.5 T at 5 K and 9 T aligned opposite to Bq. This large negative value of B-^d is consistent 
with that predicted by the resistivity measurements and supports the Jaccarino-Peter internal field- 
compensation mechanism being responsible for the origin of field-induced superconductivity. 

PACS numbers: 74.70.Kn, 76.60.-k, 75.20.Hr 



Correlations between conduction electrons and local 
magnetic moments in condensed matter physics are of 
considerable interest in situations where the properties 
of the conduction electrons are significantly tuned by the 
internal field generated by the local magnetic moments. 
These interactions can lead to a rich variety of phases, in- 
cluding superconductivity, density waves, and magnetic 
ordering. Many examples include low-dimensional or- 
ganic conductors that have been synthesized in recent 
decades 0, [^, 0] ■ It is widely accepted that their phys- 
ical properties are largely determined by the interaction 
between the donor HOMO (highest occupied molecular 
orbitals) band molecules and the anions 5, 6, 7]. 

An important example is the quasi-two dimensional 
(2D) triclinic (space group PI) salt, A-(BETS)2FeCl4, 
where BETS is bis(ethylenedithio)tetrasclenafulvalene 
(CioS4Se4H8) d, 3, S, M, Below an apphed mag- 

netic field (Bq) of about 11 T, as the temperature (T) 
is lowered it has a transition from a paramagnetic metal 
(PM) to an antiferromagnetic insulating (AFI) phase. At 
higher fields and low T there is a PM to field-induced su- 
perconducting (FISC) phase Si, Hi- 

A mechanism proposed for the FISC phase in A- 
(BETS)2FeCl4 is the Jaccarino-Peter (J-P) compensation 
effect I 111 operating in a two-dimensional (2D) system 
[^, [12, [13 1 illustrated in Fig. 1 (a) ^14]. In this model, 
the total (negative) exchange interaction (tt— d interac- 
tion, exchange constant Jjrd) between the paramagnetic 
3d Fe^"*" moments {guBSd) {g is the Lande g-factor of the 
Fe^"*" and /is is the Bohr magneton) and the conduction 
TT-electrons at the Se sites in the BETS molecule [1, [13] 
generates a large exchange field (B^^j) at the Se electrons 
aligned opposite to Bq given by 
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FIG. 1: (a) Cartoon of the interactions causing the Jaccarino- 
Peter (JP) compensation mechanism. (b) Sketch of the 
needle-like shape A-(BETS)2FeCl4 sample for the measure- 
ment. 



where < > is the average value of the Fe^+ spin polar- 
ization, and is the Lande g-factor for the 7r-electrons. 
When Bq is aligned parallel to the ac plane, the orbital 
pair breaking effect for the 7r-electrons is minimized. At 
low T (T < 5 K) for 17 T < Bq < 45 T [in the FISC 
phase] [3|, < Sd > is nearly saturated and it is expected 
that ]B^d] ~ 33 T [l2]. Thus, the magnitude of the effec- 
tive field at the Se electrons, JBq — Bird] = ]Bo] — 33 T, 
is small enough to permit the FISC phase. Also, for jBg] 
< ]B7r(i] and JBq] > ]B7rd] the spin polarization of the 
conduction electrons is respectively antiparallel and par- 
allel to Bq, a feature that can be probed by the hyperfine 
frequency shift of the ''^Se NMR signal. 

This description in terms of the J-P mechanism is sup- 
ported by the fact that its iso-structural nonmagnetic 
and non-3d-electron analog A-(BETS)2GaCl4 exhibits a 
behavior fl^, [l^ that is completely different from that of 
A-(BETS)2FeCl4. Even though the above model of FISC 
in A-(BETS)2FeCl4 is widely accepted, it needs further 
experimental confirmation. 



Nuclear magnetic resonance (NMR) is a versatile local 
probe that is capable of directly measuring the distri- 
bution of internal magnetic field and the electron spin 
dynamics on the atomic scale. Thus, it can be used as a 
tool to test the validity of the J-P mechanism for FISC 
in A-(BETS)2FeCl4. 

""^Se-NMR measurments jlTj] have been reported for 
a single crystal of A-(BETS)2FeCl4 with Bq = 14.5 T 
aligned in the ac plane (PM phase) Q. But the value 
obtained for |B7rd| is 23 T, which is ~ 30% smaller than 
the 33 T predicted by the electricial resistivity measure- 
ments 12 3] and a theoretical estimate . Also, these 
measurements do not include other alignments for Bq. 

In this paper, we report ^^Se-NMR spectrum and 
frequency shift measurements in a single crystal of A- 
(BETS)2FeCl4, for 2.5 K < T < 30 K over a range of 
alignments of Bq = 9 T in the plane ± to the c-axis. At 
this value of Bq, the PM-AFI transition is at 3.5 K. Anal- 
ysis of these results gives IBtt^I ~ (32.7 ± 1.5) T aligned 
opposite to Bq at low T and i?o > 9 T where the Fe"^"*" 
magnetization is almost fully saturated (the saturation 
value of < S'd > is slightly less than 2.5 at 9 T due to 
Fe'^"'' — Fe^+ antiferromagnetic interaction [31), consis- 
tent with the predicted value 33 T 0, [H, . It supports 
the J-P compensation as the mechanism for the FISC 
phase in A-(BETS)2FeCl4. An important input for this 
work is the Fe^^ magnetization that obtained from pro- 
ton NMR measurements on A-(BETS)2FeCl4 [Hi, which 
are not sensitive to conduction electron contributions. 

The sample used for these measurements was grown us- 
ing a standard method without ^^Se enrichment (^^Se 
natural abundance = 7.5%). Its dimensions are a* x h* 
X c = 0.09 mm x 0.04 mm x 0.80 mm [Fig. 1 (b)], cor- 
responding to a mass oi ^ 7 jig with ~ 2.0 x 10^^ ^^Se 
nuclei. Because of the small number of spins, a small 
microcoil with a filling factor (~ 0.4) was used. For most 
acquisitions, 10'*-10^ averages were used on a time scale 
of ~ 5 min for 10"* averages. The gyromagnetic ratio of 
^^Se, 7 = 8.131 MHz/T, is used for data analysis. The 
sample and coil were rotated on a goniometer (rotation 
angle 0) whose rotation axis is along the lattice c-axis 
(needle direction), which is _L Bq. Based on the crystal 
structure the direction of the Se Pz orbital is 76.4° 
from the c-axis (20| . Then, the minimum angle between 
Pz and Bq during the rotation of the goniometer is (/)min 
= 13.6°. 

Figure 2 shows the ''^Se-NMR absorption spectrum 
(x") of A-(BETS)2FeCl4, plotted as a function of the fre- 
quency shift V — vq (i/Q = 72.90 MHz) at a few T from 
30 to 2.5 K with Bq = 9 T || a', where v — has a weak 
T-dependence as discussed in more detail later. 

These spectra measure the distribution of the local 
magnetic field at the different ''^Se nuclei in the sample. 
As discussed in more detail later, the local field respon- 
sible for this distribution is dominated by the sum of 
Bq, the dipole field from the Fe'^+ spins (Bdip), and the 
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FIG. 2: (Color online) ^'^Se-NMR absorption spectrum %" as 
a function of frequency shift aX v — uo {vo = 72.90 MHz) at 
different temperatures with Bq = 9 T || a' in A-(BETS)2FeCl4. 
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FIG. 3: (Color online) T-dependence of the '^^Se-NMR 
linewidth (FWHM) A/ (solid red circles) and the modi- 
fied Brillouin function fit of the Fe''+ magnetization 
Ma{xo + x') (solid blue line) of A-(BETS)2FeCl4 at Bo = 
9 T II a' . The error bars are our estimated uncertainty. 



hyperfine field from the Se conduction electrons, whose 
polarization is strongly influenced by Bird. The spectra 
are characterized by (1) the full width at half maximum 
(FWHM) linewidth (A/) which represents the internal 
magnetic field distribution, and (2) the frequency {v) of 
the center of x" which measures the average of the hy- 
perfine field from the conduction electrons that coupled 
to the Fe^+ ions. 

Figure 3 shows A/ as a function of T for Bq = 9 T || a' . 
Also shown is a fit to the Fe^+ magnetization Md{xo + x') 
[unit: 10'^ emu/mol.Fe] provided by the modified Bril- 
louin function [l9| Bjm{,x) = Bj{xo + x'), where Bj{x) 
is the standard Brillouin function [2l|, J = S'rf = 5/2, 
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, the d-d Fe^+ exchange field B' 
Bj{xq), and Jdd = — 1.7 K is the d-d exchange 
parameter. The fit parameters for Bjm{x) are obtained 
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FIG. 4: (Color online) (a) T-dependence of the ^^Se-NMR 
resonance frequency f (solid red) for Bo = 9 T || a' in A- 
(BETS)2FeCl4. The solid line is a fit to the modified Bril- 
louin function Fe'^"'' magnetization [l^ Md{xo + x'). (b) v vs 
Mi{xo + x') above 7 K for Bq = 9 T || a' in A-(BETS)2FeCl4. 
The error bars are our estimated uncertainty. 



from proton NMR measurements fio']. In Fig. 3, A/(T) 
increases from ~ 90 kHz to ^ 200 kHz as T is lowered 
from 30 to 5 K. Also, Md{xo + x') provides a good fit. 
Since the susceptibility of the BETS molecules is small 
and nearly independent of T, Md is the main source of 
A/(r) in the PM phase. The size of A/(T) is attributed 
to a distribution of Bdip and the hyperfine field of the 
Fe^"*" across the different Se sites. Since measurements of 
the ^^Se-NMR spin-echo_decay indicate a homogeneous 
linewidth of ~ 10 kHz 19], it follows that x" is strongly 
inhomogeneously broadened. 

A similar T— dependence is also observed in t^, which 
is plotted as a function of T in Fig. 4 (a) and Mii{xq+x') 
in Fig. 4 (b). In the PM state above 7 K, a good fit 
to V (uncertainty ± 3 kHz) is obtained using 



[73.221 - 3.0158 x 10"-'A/d(a;o 



MHz. (2) 



This result is a strong indication that the 
T— dependence of v is dominated by the hyperfine 
field from the Fe'^"'' magnetization. The negative sign of 
the contribution from Ma is very important. It indicates 
that the hyperfine field from the Fe'^"'" magnetization 
is negative, i.e., opposite to Bq, as needed for the J-P 
comensation mechanism. 

A more informative result is shown in Fig. 5, where v 
— vq {vq — 72.90 MHz) is plotted as a function of both 
(j) and T. As shown in the lower right of Fig. 5, the 
^-axis is chosen to be parallel to the c-axis, Pz is the z- 
component of the BETS 7r-electron orbital moment, the 
X-axis is in the the c-pz plane, and Bq is in the xy-plane 
(_L c) rotated by 4> from the x-axis. The angle = 0° 
corresponds to Bq || the c-pz plane and the minimum 
angle between Bq and pz is ^min = 13.6°. The solid 
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FIG. 5: (Color online) '^^Se-NMR resonance frequency z/ as a 
function of angle at several T for the rotation of Bo = 9 T 
about the c-axis in A-(BETS)2FeCl4. The error bars arc our 
estimated uncertainty. 



lines are a fit to the following theoretical model based 
upon the hyperfine coupling to the BETS Se 7r-electrons 
whose polarization is affected by the exchange field from 
the Fe"^+ magnetization. 

According to NMR theory [2^, the contributions to 
the Hamiltonian (Hi) of the ^^Se nuclear spins can be 
expressed as 



Hr = H 
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where Hjz is the Zeeman Hamiltonian of the ^^Se nuclei 
in Bq, Hii is the '"''S'e— ""^Se dipole-dipole interaction 
Hamiltonian, H^^ is the direct hyperfine coupling of the 
^^Se nucleus to the BETS 7r-clcctrons generated by Bq, 
-ffjj is the indirect hyperfine coupling via the 7r-electrons 
to the 3d Fe'^"'" spins (field B^^), and H^^^ is the dipolar 
couphng to the Fe*^"*" which gives Bdip. Here, Bdip is 
calculated [2^ using the sum of the near dipole, the bulk 
demagnetization and the Lorentz contributions 2^, 2^ . 
All of these terms contribute to the static local magnetic 
field at the ''^Se sites and all but the first cause the ''^Se- 
NMR frequency shifts. 

Based upon the lattice structure and the shape and 
size of the sample, we calculated the angular dependence 
of Bdip and estimated [lO] that the contributions of Hn 
is negligible. Thus, 
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= [B„ + BA^,,{Bo,T)][l + Kc + Ks{(t)')] 

+Ks{^')B^d{Bo,T) (5) 
~ Bo(l + i^c)+Bdip(So,T) 

+if,(0')[-Bo + S,d(Bo,T)] (6) 



3 



where (j)' is the angle between Bq and the Pz direction, 
and Kc and Ks are respectively the (orbital) chemical 
shift and the (spin) Knight shift of the BETS Se tt- 
electrons. The approximation in the second line corre- 
sponds to dropping the terms Bd_i-p{Bo,T)[Kc -\- Ks{(j)')]j 
which are the product of small quantities. It can be 
shown 17, 0, 24, that for rotation of the sample 
about the c-axis with Bq in the plane _L c, Ks{4>') is 
given by 

K,{cp') = i^i,o + ifa„(</'') (7) 
X,(0) = ii:iso + /^ax[3cOs2 0COs2,?!)„,in - 1], (8) 

where Kiso and K^a{(j)) are the isotropic (independent 
of (f) and axial (anisotropic) parts of the Knight shift, 
respectively, ii'iso(ax) is a constant determined by the 
isotropic (axial) hyperfine field [^iso(ax)] produced by the 
4p^ spin polarization of the BETS Se 7r-electrons [iTJ . 
The quantities that determine K^^ are 



A, 



and 
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where xbets = 4.5 x 10""' emu/mol is the BETS tt- 
electron susceptibility [1^, = 2, Na is the Avogadro's 
number, and = + 38.6 kOe/^B reported by S. Takagi 
et al. [251. This value of is based on theoretical 



calculations < r 
Bohr radius) 



>4p = 9.28 a"^ 
27|, and ctso = 0.166 



= 0.529 A, the 

^ ^ obtained for A- 

(BETS)2GaCl4, which has essentially the same BETS- 
molecules as A-(BETS)2FcCl4. By using these values, 
one obtains A'ax = 15.3 x 10^"*. Its uncertainty is not 
known to us and is not included in our analysis. 

The angular dependence of the shift in Ks{(j)) [Eqs. 
(7)-(8)] has been used for the fit of the 5 K data shown 
in Fig. 5 using the relation v — vq — ai[3 cos^ cos^ 13.6° 
— 1] + bi, with the fit values ai — ~ 313 kHz and bi — 
+ 221 kHz (uncertainty ± 8 kHz). 

From Eqs. (6) — (8), the formula for B^^ at a fixed Bq 
and Tq from the difference in v (Ai^) at the two angles 
and (j)2 is 

n ^^ Aiy{To,cj)l,(t)2) Bdip(0l) - -Bdip(02) 



AK: 



(11) 



where 



33 T obtained from the electrical resistivity measurement 
n , Tsj and the theoretical estimate 0] . Also, a small in- 
crease in < Srf > is expected as Bq is increased from 9 
T to 33 T at 5 K. From the modified Brillouin function 
analysis used for the proton NMR linewidth we ex- 
pect this increase in < > to be a factor 1.05 ± 0.05, 
which corresponds to an adjustment to 5^^(33 T, 5 K) 
= - 34.3 ± 2.4 T. 

Similar results for the exchange field B^rd can also be 
obtained from the T-dependence of at a fixed cj) with 
the data in Figs. 4 and 5 using the same kind of fit 
as Eq. (2). But the uncertainty in the value obtained 
for _B^d with this type of analysis is large enough that 
we do not present it here. An important test we plan 
to do in the future is to extend the measurements to 
Bq > |i3^d| to find the value of Bq = |-B,rd|, where the 
angular dependence in Fig. 5 disappears [Eqs. (6) — (8)] 
and above which the sign of Ks{(j)') changes from negative 
to positive. 

Figure 5 shows two values of 4> {4>tii and 0ti2) where 
u becomes independent of T . The measured difference 
'/'Til — 0TI2 = 65.0° ± 2.0°, which should be symmetric 
around 90°, or (/)tii = 57.5° ± 1.0° and 0ti2 = 122.5° ± 
1.0°. By using Eqs. (6) and (8) and neglecting the very 
small contribution from _Bdip (^dip << I^Trdl), it can be 
shown that the condition for (p^n (i = 1 , 2) is -f^iso + 
/Taxpcos^ (/)Tii cos^ (/imin — 1] = 0. From this relation, 
the measured values of 0tii and (^ti2, and i^ax = 15.3 x 
10-"*, one obtains K-.^o = (2.8 ± 0.7) x 10"'', or Ki^o/K^^ 
= 1/(5.5 ± 1.4). 

In summary, our results of "^^Se-NMR spectrum and 
frequency shift measurements in A— (BETS)2FeCl4 indi- 
cate that the Fe'^"'" spins have a strong antifcrromagnetic 
coupling to the BETS 7r-electrons, and wc determined 
the TT — d exchange field to be -Bjrd = •/•n-d < Sd > /qhb 
= - 32.7 ± 1.5 T at Bo = 9 T and 5 K, with an ex- 
pected value of - 34.3 ± 2.4 T at Bq ~ 33 T and 5 K. 
This large negative value of By^d (or J^^) is consistent 
with that predicted by the resistivity measurements, and 
supports the Jaccarino-Peter internal field-compensation 
mechanism being responsible for the origin of the FISC 
phase in A-(BETS)2FeCl4. 
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The conditions T = 5 K, 0i = 90°, and 02 = 0° give 
Az/(5K, 90°, 0°) = 880 ± 26 kHz (from Fig. 5) and Aifan 
= 4.42 X lO^'^ (estimated error ± 3%). Our calculated 
value of Bdip(0°) - Bdip(90°) at 5 K is 3.63x10"^ T and 
the field used was Bq — 9.0006 T. These values then give 
B^rd = — 32.7 ± 1.5 T at 5 K (ahgned opposite to Bq = 
9.0006 T), which is very close to the expected value of — 
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